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Material specific electronic band structure of the electron-doped high-Tc cuprate 
Ndi.85Ceo.i5Cu04 (NCCO) is calculated within the pseudo gap regime, using the recently 
developed generalized LDA+DMFT+Ek scheme. LDA/DFT (density functional theory within 
local density approximation) provides model parameters (hopping integral values, local Coulomb 
interaction strength) for the one-band Hubbard model, which is solved by DMFT (dynamical 
mean-field theory). To take into account pseudogap fluctuations LDA+DMFT is supplied with 
"external" k-dependent self-energy Ek, which describes interaction of correlated conducting 
electrons with non-local Heisenberg-like antiferromagnetic (AFM) spin fluctuations responsible for 
pseudo gap formation. Within this LDA+DMFT+Ek approach we demonstrate the formation 
of pronounced "hot-spots" on the Fermi surface (FS) map in NCCO, opposite to our recent 
calculations for Bi2Sr2CaCu208-i (Bi2212), which have produced rather extended region of FS 
"destruction". There are several physical reasons for this fact: (i) the "hot-spots" in NCCO are 
located closer to Brillouin zone center; (ii) correlation length of AFM fluctuations ^ is larger for 
NCCO; (iii) pseudogap potential A is stronger, than in Bi2212. Comparison of our theoretical 
data with recent bulk sensitive high-energy angle-resolved photoemission (ARPES) data for NCCO 
provides good semiquantitative agreement. Based on that comparison alternative explanation of 
the van-Hove singularity at -0.3 eV is proposed. Optical conductivity both for Bi2212 and NCCO is 
also calculated within LDA+DMFT+Ek and compared with experimental results, demonstrating 
satisfactory agreement. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.27+a, 71.30.+h, 74.72.-h 



I. INTRODUCTION 

There is a good reason to think that proper descrip- 
tion of pseudogap regime is the avenue of approach to 
the physical nature of high-Tc superconductivity^ . An- 
gle resolved photoemission spectroscopy (ARPES) has 
been coming along this way very well in recent years. 
One of the test compounds for ARPES is the hole doped 
Bi2Sr2CaCu208-5 (Bi2212) system. Another example is 
the electron doped high-T^ Nd2-:rCe:rCu04 (NCCO). To 
this end there are plenty of experimental ARPES data on 
Bi2212 and NCCO (see the review^). Fermi surface (FS) 
maps, quasiparticle band dispersions and even self-energy 
lineshapes within mapping on some models are reliably 
extracted from modern ARPES data2.. 

There are several interesting physical phenomena asso- 
ciated with pseudo gap regime (in the normal underdoped 
phase) : partial FS "destruction" and folding of band dis- 
persions (shadow bands) for both compounds Bi2212 and 
NCCO^. Despite evident similarities of experimental ob- 
servations for these two systems there is one striking dis- 
tinction. The FS of Bi2212 has so called Fermi "arcs" 
around (7r/2,7r/2) point (looking like a part of noninter- 
acting FS), while towards BZ borders sharply defined FS 
just vanishes. In its turn NCCO also possesses slightly 
degraded Fermi "arcs" , but in the vicinity of BZ borders 
noninteracting FS is almost restored. In between there 
are well known "hot-spots" — areas of FS "destruction" 



around the points where FS intersects umklapp BZ bor- 
der. These "hot-spots" are not observed so obviously for 
Bi2212. The aim of the present paper is to show the 
origin of this NCCO "hot-spots" behavior. 

Both systems under consideration are usually treated 
as Mott insulators at moderate doping or, in other words, 
as strongly correlated metals. Modern technique to 
solve the Hubbard model is dynamical mean-field the- 
ory (DMFT), which is exact in infinite dimensions^. 
However, the quasi two-dimensional nature of high-Tc 
compounds is well known. To overcome the local na- 
ture of the DMFT approximation, we have proposed 
recently a semiphenomenological DMFT+Sk computa- 
tional scheme^i^i^. Here an additional self-energy Sk de- 
scribes non-local correlations induced by (quasi) static 
short-ranged collective Heisenberg-like antiferromagnetic 
(AFM) spin (pseudogap) fluctuations'^^:-. Assuming addi- 
tive form of the self-energy within DMFT+Sk approach 
one can preserve conventional DMFT self-consistent set 
of equations. To take into account material specific prop- 
erties of both Bi2212 and NCCO we perform first prin- 
ciple one-electron density functional theory calculations 
within local densty approximation (DFT/LDA)^. Then 
LDA results are incorporated into DMFT+Sk in accord 
with LDA+DMFT ideologylA. To solve the effective 
single impurity problem of the DMFT we employ here 
the reliable numerical renormalization group approach 
(NRG)iiii^. Such combined LDA+DMFT+Sk scheme 
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by construction is particularly suitable to describe elec- 
tronic properties of real high-Tc materials at finite doping 
in the normal state. 

The DMFT+Sk approach was extensively applied dur- 
ing the last years to describe formation of pseudogap 
within strongly correlated metallic regime of single-band 
Hubbard model on the square lattice^i^i^. We have also 
generalized DMFT-|-Ek approach for the account of static 
disorder effects^. Later we have derived DMFT-f Ek ap- 
proach to calculate two-particle properties (such as op- 
tical conductivity)^'^. Recently DMFT-|-Sk was used by 
us to analyze the general problem of metal - insulator 
transition in strongly disordered and strongly correlated 
systemsii. 

LDA-I-DMFT-I-Sk scheme was already used to describe 
pseudogap regime in "realistic" calculations for Bi2212i^. 
In the present paper we apply this approach to NCCO 
aiming to describe the characteristic differences of its 
electronic structure as compared to Bi2212. 

The paper is organized as follows. In section |TT] 
we present a short introduction into an ab initio 
self-consistent generalized combined LDA+DMFT-I-Sk 
scheme and its extension for two particle properties (op- 
tical conductivity). Section IIIII contains Bi2212 and 
NCCO material specific information: LDA calculated 
band structure, Fermi surfaces and details on some 
model parameters calculations. Results and discussion 
of LDA-hDMFT+Sk calculations for Bi2212 and NCCO 
and comparisons with experimental data are presented 
in section IIVI Section |V] concludes our manuscript with 
the summary and discussion of some of the remaining 
problems. 



II. COMPUTATIONAL METHOD 

To introduce spatial length scale (nonlocal correla- 
tions) into conventional DMFT method'^ we have recently 
proposed the generalized DMFT+Sk approach^i^, with 
computational scheme, shown in Fig. [1] which contains 
the flow diagram of self consistent DMFT-I-Sk set of 
equations. First we guess some initial local (DMFT) 
electron self-energy Ti{iuj). Second we compute by any 
available technique (for the chosen model) k-dependent 
"external" self-energy Ek(«a;) which can be in general 
case a functional of S(ja;). Then neglecting interference 
effects between the self-energies (which in fact is the ma- 
jor assumption of our approach) we can set up and solve 
the lattice problem of DMFT (step 3 in Fig.[T]). Then at 
step 4 we define effective Anderson single impurity prob- 
lem which is to be solved by any "impurity solver" to 
close DMFT-I-Ek equations. 

The additive form of self-energy (at step 3 in Fig. [T|) is 
in fact an advantage of our DMFT-HEk approachii^i^. It 
allows one to preserve the set of self-consistent equations 
of standart DMFT'^. However there are two distinctions 
from conventional DMFT. During each DMFT iteration 
we recalculate corresponding k-dependent self-energy 
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Solution of effective single impurity Anderson problem 



FIG. 1: Flow diagram of DMFT-fEk selfconsistent loop, ii 
corresponds to lattice problem and d - impurity problem vari- 
ables. 



I]k(A*,^, [^('^)]) within some (approximate) scheme, e.g. 
taking into account interactions with collective modes or 
order parameter fluctuations, and the local Green func- 
tion Gu{iuj) is "dressed" by Ek at each step. When input 
and output Green's functions (or self-energies) converge 
to each other (with prescribed accuracy) we consider the 
obtained solution to be self-consistent. Physically it cor- 
responds to the account of some "external" (e.g. pseudo- 
gap) fluctuations, characterized by an important length 
scale ^, into fermionic "bath" surrounding the effective 
Anderson impurity of the usual DMFT. 

In the present work Ek(i^) represents interaction of 
correlated electron with antiferromagnetic (AFM) pseu- 
dogap fluctuations. To calculate Ek(w) for the case 
of random fleld of pseudogap fluctuations (assumed to 
be (quasi) static and Gaussian, which is valid at high 
enough temperatures^i^) with the dominant scattering 
momentum transfers of the order of characteristic vector 
Q = (7r/a,7r/a) (a - lattice parameter), typical for AFM 
fluctuations ("hot-spots" model^), the recursion proce- 
dure proposed in Refsi^'^'^ is used, with material speciflc 
generalizations as described in detail in Refs.^^. 

There are two important parameters characterizing 
pseudogap regime in our scheme: the pseudogap energy 
scale (amplitude) A and spatial correlation length S^^. 
Actually we prefer to take A and ^ determined some- 
how from experiments. However, we can also use cer- 
tain model estimates to calculate them microscopicallj«^. 
Both approaches are used below. 

To calculate optical conductivity we use our gener- 
alization of DMFT+Ek for calculation of two - parti- 
cle properties (vertex parts) as described in detail in 
Ref.^^, with material dependent parameters provided by 
LDA-fDMFT-|-Ek and vertex corrections due to pseudo- 
gap fluctuations calculated using the recursion relations 
derived earlier in Ref 
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FIG. 2: LDA bands (gray lines) for Bi2212 (left) and NCCO 
(right) along BZ high symmetry directions. For both pan- 
els solid black lines correspond to effective x^ — y^ symmetry 
Wannier-like state dispersions. Zero energy corresponds to 
Fermi level. 



III. LDA BANDS AND FS OF NCCO AND 
BI2212, EFFECTIVE MODEL PARAMETERS 

As a first step of our LDA+DMFT+Ek hybrid scheme 
we perform LDA band structure calculations. For both 
compounds ideal tetragonal bcc crystal lattice with space 
symmetry group I4/mmm is reported (for Bi2212 see 
Ref,il, for NCCO RefJ^). Physically relevant structural 
motif for high-Tc materials is Cu02 plane. There are two 
Cu02 planes displaced close to each other in the unit cell 
of Bi2212, while NCCO has just one. We have done LDA 
calculations of electronic band structure within the lin- 
earized muffin-tin or bital (LMTO) basis set?0. Rcsuhs 
are presented as gray lines in Fig. [21 left panel — Bi2212 
bands, right — NCCO. Our band structures agree well 
with previous works Ref.^ and Ref.-- for Bi and Nd com- 
pounds correspondingly. 

To calculate hopping integral values for Bi system 
Wannier functions projecting method^'^ in the LMTO 
framework^^ was applied. Hopping integrals of Nd com- 
pound were obtained by using the so called NMTO 
method^^ (see Table IIIII) . Values of hopping integrals 
computed by these two methods agree well for the respec- 
tive compounds^. In Fig. [5] by black line represents the 
dispersion of effective — Wannier-like orbital which 
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FIG. 3: LDA calculated Fermi surfaces for Bi2212 (left) and 
NCCO (right) in the quarter of the Brillouin zone. Diagonal 
line corresponds to (AFM) umklapp scattering surface. 

crosses the Fermi level and is most interesting physically. 
These dispersions correspond to hopping integral values 
(tight-binding parameters) given in Table IIIII 

Fig. [3] contains noninteracting LDA Fermi surfaces 
(FS) in {kx,ky) plane for quarter of the first Brillouin 
zone (BZ) (left panel - Bi, right panel - Nd systems). 
The shape of these FS is defined by the tight-binding 
parameters from Table Hill Diagonal line corresponds to 
the AFM folded BZ border. In the left panel for Bi2212 
one can see two FS sheets. It is caused by finite hopping 
between two neighboring Cu02 layers — the so called 
bilayer splitting. The value t± is listed in the Table Hill 
Simplest tight-binding expression for the bilayer splitting 
derived by Andersen et al^ is used in our calculations 
(for details see RefJ^). 

It is important to note the "hot-spots" positions (in- 
tersections of FS with AFM umklapp surface) for both 
materials. For Bi2212 it is (0.47,2.66)7r/a and for NCCO 
(0.95,2. 19)7r/a. Whereof one can see that for NCCO the 
"hot-spots" are located farther away from the BZ border 
than in Bi2212. Let us remind that pseudogap fluctua- 
tions scatter electrons from the vicinity of one "hot-spot" 
to the vicinity of another "hot-spot" , i.e. by the scatter- 
ing vector of the order of Q. The effective scattering area 
around the "hot-spot" is defined by inverse correlation 
length of these fluctuations. Let us also remember 
that (n/afi) point is surrounded by four BZ from dif- 
ferent sides. Consequently if the "hot-spot" is closer to 
(7r/a,0) point and ^ is small enough, FS is "destroyed" in 
rather wide region close to BZ border crossings. Thus one 
can expect for NCCO "hot-spots" to be observed more 
explicitly, while the part of the FS close to the BZ bor- 
der crossings be less affected by pseudogap fluctuations 
in contrast to Bi2212. 

The values of local Coulomb interaction U for — 
orbital were obtained by constrained LDA method'^* (for 
values see Table IIIII) . These values are of the order of 



TABLE L Calculated energetic model parameters for Bi2212 
and NCCO (eV). First four Cu-Cu in plain hopping integrals 
t, t' , t" , t'" , interplain hopping value t±, local Coulomb in- 
teraction U and pseudogap potential A. 





t 


t' 


t" 


t'" 


t± 


U 


A 


Bi2212 


-0.627 


0.133 


0.061 


-0.015 


0.083 


1.51 


0.21 


NCCO 


-0.44 


0.153 


0.063 


-0.0096 




1.1 


0.36 
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FIG. 4: LDA+DMFT+Ek spectral functions contour plot for 
Bi2212 (upper panel) and NCCO (lower panel) along BZ high- 
symmetry directions of first Brillouin zone r(0, 0) — X(7r, 0) — 
M(7r, tt) — r(0, 0). Zero energy corresponds to Fermi level. 
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FIG. 5: Comparison of effective LDA — bands (upper 
panel) and LDA+DMFT+Ek quasiparticle dispersions (lower 
panel) Bi2212 and NCCO along BZ high-symmetry directions. 
Zero energy corresponds to Fermi level. 



2-3 t for both systems. It is in fact quite a bit smaller 
than the values many people believe should be used in a 
model calculations (usually about 4-6 t, see e.g. Refi^). 
At the same time, our previous experience with con- 
strained LDA computation shows that it gives reason- 
able estimates for the U value in a number of other ox- 
ides^. However, to analyze further the influence of U 
value on observable quantities we performed additional 
LDA+DMFT+Sk computations for increased values of 
U . The short discussion of the results will be given in 
Sec. El Values of A for both systems were calculated 
as proposed in Ref. — (Table HII)) . (for more details see 
Appendix |^ . Correlation length ^ was taken from ex- 
periments, i.e. ^ ~ 10a for Bi2212^ and ^ ^ 50a for 
NCCO^. 



IV. NCCO VS. BI2212 LDA+DMFT+Ek 
RESULTS AND EXPERIMENTAL DATA 

A. Quasiparticle dispersions 

Finite temperatures and interactions lead in general 
to finite life-time effects. Thus, instead of quasiparticle 
dispersions expressed by the usual dispersion curves (like, 
for example, in DFT/LDA) in Fig. 2] we display contour 
plots of corresponding spectral functions A(a;,k): 

A(w,k) = --ImG'(w,k), (1) 

where G(a;,k) is the retarded Green's function obtained 
via our LDA+DMFT+Ek scheme (shown in Fig. [J) with 
appropriate analytic continuation to real frequencies. 

Prima facie both compounds Bi2212 (left panel of 
Fig. H]) and NCCO (right panel of Fig. 0]) have similar 



quasiparticle bands. There are two bands in each case 
instead of just one in case of DFT/LDA. Of these, the 
most broad and intensive band predominantly follows the 
noninteracting DFT/LDA band (see Fig. [2]). The second 
band in our case is an AFM-like reflex (shadow band) of 
the latter one and has much weaker intensity. This is a di- 
rect effect of self-energy Ek due to pseudogap fluctuations 
introduced into the conventional DMFT (for discussion 
see Ref,^). 

As we discussed above, finite life-time (interaction) ef- 
fects should be especially strong around the (7r/a,0) point 
(AT-point in Fig. 2]) . This is clearly visible in both panels 
of Fig. m However in Sec. [m] we showed that Bi2212 
"hot-spot" is much closer to the AT-point. Thus we see 
large quasiparticle band broadening (do not forget bi- 
layer splitting effects) on the Fermi level. For NCCO 
there is no sizeable broadening on the Fermi level close 
to X-point. Both branches go below Fermi level at about 
-0.5 eV. In NCCO the Fermi "arc" is much closer to umk- 
lapp surface, so pseudogap effects are significant around 
(7r/2a,7r/2a) point (middle point of the V — M direction) 
in contrast to Bi2212. 

Fig. [5] shows changes from effective —y^ LDA bands 
(left side) to LDA+DMFT+Sk quasipartical bands 
(right side) for both Bi2212 and NCCO. Quasiparticle 
bands on the right panel of Fig. [S] represent the positions 
of the maxima of spectral functions shown in Fig. [D In 
Bi2212 the shadow band and the quasiparticle band in- 
tersect each other in the "hot-spot" close to the X-point. 
While in NCCO there is no such intersection but the 
shadow and quasiparticle bands are quite parallel around 
AT-point. Close to the (7r/2a,7r/2a) point we observe a 
kind of precursor of the dielectric AFM gap. Nothing of 
that sort is observed for Bi2212. Also one should men- 
tion that calculated shadow band in Bi2212 is actually 
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FIG. 6: LDA+DMFT+Ek spectral functions for Bi2212 (up- 
per panel) and NCCO (lower panel) along of noninteracting 
FS in 1/8 of BZ. Dashed- black line corresponds to "hot-spot" . 



an order of magnitude less intensive than in NCCO. 



B. Spectral functions 

Fig. M displays LDA+DMFT+Sk spectral functions 
((!]) along 1/8 of noninteracting FS from the nodal point 
(top curve) to the antinodal one (bottom curve) (A and 
B points in Fig. [31 correspondingly). Data for Bi2212 is 
given in left panel, NCCO — right panel of Fig. \6\ For 
both compounds antinodal quasiparticles are well-defined 
— sharp peak close to the Fermi level. Going to the 
nodal point quasiparticle damping grows and peak shifts 
to higher binding energies. This behavior is confirmed 
by experiments Refsi^i^i (for brief comparison with ex- 
periment see Refj^). Again there are some differences 
between these two compounds. As we said before "hot- 
spots" for NCCO are closer to the BZ center. In Fig. [5] 
one can see it from the position of the dashed-black line 
which corresponds to the "hot-spots" k-point. Thus an- 
other explanation of the peaks can be given. Namely, 
for Bi2212 nodal quasiparticles are formed by low energy 
edge of pseudogap. While for NCCO they are formed by 
higher energy pseudogap edge. Also in NCCO there is 
obviously no bilayer splitting effects seen for Bi2212 (left 
panel of Fig. 

C. Comparison with ARPES data 

In Fig. El LDA+DMFT-I-Sk FS maps on the quarter 
of BZ for Bi2212 (upper left) and NCCO (upper right) 
are presented. The upper parts of Fig. [7| are just a con- 
tour plot of the spectral functions from Fig. [SI Above 
mentioned significant FS "destruction" because of pseu- 
dogap fluctuations close to the borders of BZ is clearly 
seen for Bi2212. Quite the contrary NCCO FS is almost 
restored in the vicinities of the BZ border. Vice versa 
Fermi "arc" of Bi2212 is quite sharp while for NCCO is 
rather degraded. That is afresh the consequence of the 




FIG. 7: LDA-I-DMFT-I-Ek Fermi surfaces for Bi2212 (upper 
left panel) and NCCO (upper right panel) 1/4 of BZ {kx, ky in 
units of n/a). Experimental FS for Bi2212 (lower left panel, 
ReiM) and NCCO (lower right panel, Ref.^). 



"hot-spots" to be closer to BZ center for NCCO. Slightly 
bigger value of pseudogap potential A also works towads 
Fermi "arc" smearing in NCCO. It is significant to say 
that shadow FS are come to hand. Shadow FS is found 
to be more intensive for NCCO. 

Qualitatively speaking very similar FS shapes are ob- 
served experimentally for both Bi Refi^ and Nd Refi^ 
compounds (lower parts of Fig. [7]). To our opinion such 
FS maps have material specific origin. Namely LDA cal- 
culated FS of NCCO has more curvature (left panel of 
Fig. [3]) and intersects the BZ boundary away from the 
(jr/afi) thus remaining nearly noninteracting one. While 
Bi2212 FS comes to BZ border much closer to the {ir/afi) 
point (right panel of Fig. [3]). Because of that in Bi2212 
"hot-spots" are not seen. They are spread by strong 
pseudogap scattering processes near (n/afl) point. Also 
bigger correlation length for NCCO is favorable for more 
evident "hot-spots" . 

In Fig.[5|we present LDA+DMFT-|-Ek data in compar- 
ison with recent high-energy bulk sensitive angle-resolved 
photoemission data of Ndi.85Ceo.i5Cu04. For the details 
of experiment see Ref.-"^. Lower panel of the Fig.[8]shows 
intensity plots along the high symmetry lines for NCCO 
obtained by the high-h^ ARPES. Upper panel of Fig. [Hj 
is part of the Fig. [D Here we should say that to get bet- 
ter agreement with experiment theoretical Fermi level is 
changed by 0.2 eV. 

Indeed we see quite a good agreement of 
LDA-fDMFT-l-Ek and experimental data. For the 
M — r direction there is not very much going on. 
Basically we see both in theory and experiment very 
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FIG. 8: Comparison of LDA+DMFT+Sk spectral functions 
(upper panel) for NCCO along BZ high-symmetry directions 
with experimental ARPES (Ref.^^) (lower panel). 



intensive quasiparticle band. For the M — T direction low 
intensive shadow band is not resolved in the experiment. 

More interesting situation is observed for T — X — M 
directions. At F-point there is band in the experiment 
starting at about -1.2 eV. It is rather intensive and goes 
up in energy. Suddenly there is almost zero intensity at 
about -0.3 eV. Then in the vicinity of the X-point in- 
tensity rise up again. In the X — M direction around 
-0.3 eV on the right side of X-point there is also quite 
intensive region. At a first glance one can think that it 
is the same band with matrix element effects governing 
intensity. However looking at the right panel Fig. O (see 
corresponding discussion in Sec. IIV Ap one can realize 
that this low intensity region is the forbidden gap be- 
tween shadow and quasiparticle bands. The "horseshoe" 
around AT-point is formed by the shadow band on the left 
and the quasiparticle band on the right for upper branch 
and other way round for the lower branch. It is also easy 
to see this in Fig. 2] and the upper panel of Fig. [51 As 
a consequence of that there is also intensive shadow FS 
sheets around (7r/a,0) point on Fig. [7] (upper right panel). 
Rather intensive nondispersing states at about -1.0 eV 
within experimental data can be presumably associated 
with the lower Hubbard band and possible admixture of 
some oxygen states. Let us also suppose that high in- 
tensity at -0.3 eV for X point may be interpreted not as 
a van-Hove singularity of bare dispersion^- but rather of 
high-energy pseudogap branch. 



D. Comparison with optical data 

Our recent generalization of LDA+DMFT+Sk scheme 
with respect to two-particle properties''^ allows us to an- 
alyze optical conductivity of the Bi and Nd materials 
under consideration. In Fig. [9] comparison of experimen- 
tal data with part of optical conductivities for NCCO 
(left panel) and Bi2212 (right panel) is presented. The 
way of computation is described above. Here we can re- 
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FIG. 9: Comparison of LDA-|-DMFT-|-Ek calculated optical 
conductivity spectra for NCCO (left panel) with experimental 
data (circles) Ref.^®. Solid line — theoretical results for cal- 
culated pseudogap value A=0.36 eV (dashed line corresponds 
to experimental A=0.2 eV). On the right panel there is the 
same quantity but for Bi2212 and experiment of Ref.— . 



port qualitative agreement of our theoretical curve for 
NCCO with calculated A=0.36 eV (solid line) with ex- 
periment2S. Nevertheless we find calculated pseudogap 
value to be overestimated. To improve the agreement 
we also calculated optical conductivity for experimental 
value of A=0.2 eV— (Fig. [HI dashed hue). Concerning 
Bi2212 optical conductivity (Fig. [51 right panel) one can 
point out that there is no particular structure neither 
in theory nor in the experimental data^. For Bi2212 
agreement between experimental and theoretical curves 
is reasonable. Let us note that according to our calcu- 
lations of the quasiparticle bands, spectral functions and 
FS maps are not strongly modified in case of A=0.2 eV. 



E. Influence of U value on DMFT+Sk results 

Before we summarize our results, we would like to 
discuss how the Hubbard interaction U value affects 
DMFT-fEk results, namely, the observable physical 
quantities. This question arises from ongoing scientific 
discussion in the literature. Note that our constrained 
LDA calculated U value is of the order of 2-3t. At the 
same time, it is commonly believed that U value should 
be of the order of 4-6*22ii^. To this end we have per- 
formed additional DMFT+Sk computations, using these 
values of U. 

If we take just higher values of U within DMFT+Sk 
approach, without any change of other parameters of our 
model, we obtain stronger uniform quasiparticle damp- 
ing. Spectral functions become slightly more blurred and 
FS less sharply defined. Also with the increase of the U 
value quasiparticle mass grows up a little bit. In optical 
conductivity the pseudogap anomaly also appears to be 
more damped, as we mentioned earlier in Ref.—. At the 
same time the general agreement with experiments stays 
reasonable. This is not very surprising since these values 
still belong to small or moderate correlations {U is less 
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than band width W = 8t). 

But actually, we are trying to express the pseudogap 
potential A via U (see Appendix [A|). If this connection 
is taken into account FS maps do not differ very much 
from those obtained above. However, the pseudogap ef- 
fects are getting stronger both around "hot-spots" and in 
the vicinity of (7r/2a, n/2a) k-point. In this case compar- 
ison with ARPES data for quasiparticle bands becomes 
much more worse. Especially, the bigger U values spoil 
the accord with optical data. Thus we can conclude, 
that constrained LDA calculated value of U (together 
with corresponding value of A) allows one to describe 
ARPES experimental data reasonably well, though in 
optical conductivity the size of the pseudogap is a bit 
overestimated. In other respects, we do not observe any 
qualitative changes in our results as U increase from 3t 
to U or 6t. 



CONCLUSION 



hospitality. 
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FIG. 10: Filling dependence of the pseudo gap potential A 
calculated with DMFT(QMC) and DMFT(NRG) for varying 
Coulomb interaction ([/) at temperature T=QAt on a two- 
dimensional square lattice with i'/i=-0.4. 



To summarize, the origin of evident "hot-spots" on 
NCCO FS in pseudogap regime is attributed to the de- 
tails of noninteracting electronic band structure of this 
compound. All differences in physical properties cal- 
culated within the LDA-f DMFT-I-Sk approach (quasi- 
particle bands, FS maps and ARPES data) are deter- 
mined by the fact that in NCCO "hot-spots" lie closer 
to BZ center than in Bi2212. Also stronger AFM long 
range ordering tendency in NCCO favors the clearly vis- 
ible "hot-spots" . Besides that, the qualitative behavior 
of both electron doped NCCO and hole doped Bi2212 
high-Tc systems is almost the same. We have also in- 
terpreted the new ARPES experimental data for NCCO 
by LDA-|-DMFT-|-Sk calculated quasiparticle bands and 
proposed the new mechanism for the origin of van-Hove 
like singularity at -0.3 eV. 

The results obtained further support the 
LDA+DMFT-I-Sk to be an effective method to in- 
vestigate the electronic structure of strongly-correlated 
systems. 
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APPENDIX A: PARTICLE-HOLE ASYMMETRY 
OF THE PSEUDOGAP POTENTIAL A 

Pseudogap energy scale (amplitude) A was calculated 
in Ref.5 via DMFT(QMC - Quantum Monte Carlo) sim- 
ulations for the hole doped region. Here we present sim- 
ilar results for the electron doped case obtained with 
DMFT(NRG - Numerical Renormalization Group). Us- 
ing the two-particle selfconsistent approach of Ref.— , 
with the approximations introduced in Refs.^'-, for stan- 
dard Hubbard model one can derive following micro- 
scopic expression for A: 

where we take into account only scattering by antifer- 
romagnetic spin fluctuations. The different local quan- 
tities here, such as density n and double occupancy 
< rii^nii > etc., can easily be calculated within the stan- 
dard DMFT"^. As impurity solver we have used NRG. 

In Fig. [TO] we show our results for A for both electron 
and hole dopings. One can immediately see a remarkable 
(up to one order of magnitude) particle-hole asymmetry 
of the Eq. (|Aip for large values of U. For values of U 
below or equal to 8t (which corresponds to week or mod- 
erate coupling case) this A particle-hole asymmetry is 
about factor of two. Basically in case of the Eq. (|A1[) 
it comes from particle hole asymmetry of double occu- 
pancy < niijUii > value which is enhanced as U value 
increases. These results agree with experimental obser- 
vation for pseudogap effects to be stronger for electron 
doped systems. 
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